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Corrections of orders and are calculated in the fine structure interval AE^^ = 
E{2P^I2) — E{2Pi/2) and in the hyperfine structure of the energy levels 2Pi/2 and 
2P3/2 in muonic hydrogen. The obtained numerical values AE^^ = 8352.08 fieV, 
A&f'{2Pii2) = 7964.36 ^eF, A^^/^(2P3/2) = 3392.59 ^leV can be considered 
as a reliable estimate for the comparison with corresponding experimental data and 
for the extraction of the experimental value of the Lamb shift {2P — 23) in muonic 
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I. INTRODUCTION 

Simple atoms play important role in the check of quantum electrodynamics (QED), the 
bound state theory and precise determination of fundamental physical constants (the fine 
structure constant, the lepton and proton masses, the Rydberg constant, the proton charge 
radius, etc) [l|, On the one hand the essential progress achieved in the last years in a 
more precise determination of fundamental physical parameters |13( is closely related with 
the growth of the theoretical accuracy in the calculation of the fine and hyperfine structure 
of the energy levels of hydrogenic atoms, the lepton anomalous magnetic moments (AMM) 
and on the other hand with the perfection of experimental methods in the investigation of 
atomic spectra of the hydrogen, muonium, positronium Light muonic atoms [1,H,0| 

(muonic hydrogen (/xp), muonic deuterium, ions of muonic helium etc.) are distinguished 
among simple atoms because of the strong influence of the vacuum polarization (VP) effects 
H, Ipl , the recoil effects [l^ , the nuclear structure and polarizability effects [l^, EH 
EB; HH ^3 the structure of their energy levels. So, the comparison of the theoretical 
value of the Lamb shift (2P — 25*) in muonic hydrogen with experimental data from the 
Paul Sherrer Institute (PSI) [1^ E^ will make it possible to obtain a more precise value 
of the proton charge radius and to check the quantum electrodynamics with the accuracy 
10~^. Considering that the energy interval {2^P^i2 — 2^5*1/2) is investigated at PSI with an 
accuracy 30 ppm, it is important to calculate precisely not only the Lamb shift value, but 
also the fine and hyperfine structure of the S - wave and P - wave energy levels in the atom 

The ene rgy levels of light muonic atoms were theoretically studied many years ago in 
[20I l2l[ 2^ . 231 2^ 1 both on the basis of the Dirac relativistic equation and the nonrelativistic 
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approach, taking into account corrections by the perturbation theory (PT). In these papers 
(see also references jisl . |2^) the basic contributions to the fine and hyperfine structure of 
the P- energy levels in muonic hydrogen were evaluated with the accuracy 0.001 meV . The 
structure of 5'-wave states in (/ip) was studied in ji^, IH, |2^ 3^1 accounting corrections of 



order and a^. In this work we continue the investigation of the energy spectrum of (/xp) 
in the P-wave part. The aim of the present study is to calculate such contributions of order 
and a® both in the fine and hyperfine structure of the energy states 2Pi/2, 2P3/2, which 
are connected with the electron vacuum polarization, the recoil effects, the muon anomalous 
magnetic moment and the relativistic corrections. The role of all these effects is crucial 
to attain the desirable accuracy. Our goal also consists in the refinement of the earlier 
performed calculations in 2l|, |23| and in the derivation of the reliable numerical estimate for 



the structure of P-wave levels in the atom (/ip), which can be used for the comparison with 
experimental data. Modern numerical values of fundamental physical constants are taken 
from Ref.j^: the electron mass rrie = 0.510998918(44) ■ 10^^ GeV, the muon mass = 
0.1056583692(94) GeV, the fine structure constant a'^ = 137.03599911(46), the proton mass 
rrip = 0.938272029(80) GeV, the proton anomalous magnetic moment k = 1.792847351(28), 
the muon anomalous magnetic moment = 1.16591981(62) ■ 10~^. 



II. FINE STRUCTURE OF P - WAVE STATES 



Our approach to the investigation of the energy spectrum of muonic hydrogen is based 
on the use of quasipotential method in quantum electrodynamics [31I. l3^ , where the two- 
particle bound state is described by the Schroedinger equation. The basic contribution to 
the muon and proton interaction operator is determined by the Breit Hamiltonian [33]: 
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where mi, m2 are the muon and proton masses, fi = mim2/ {mi + 7712) is the reduced mass, 
AV^^ is the muon spin-orbit interaction: 
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/^yhfs jg ^YiQ proton spin-orbit interaction and the interaction of the muon and proton spins. 
The leading order {Za^ contribution to the fine structure is determined by the operator 
AV^^. As it follows from the Eq.(2), AV^^ includes also the recoil effects (the Barker-Glover 
correction fl^) and the muon anomalous magnetic moment correction. Let us write the 
fine structure interval for the atom (fip) in the form: 
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FIG. 1: The contribution of order a(Za)^ to the fine and hyperfine structure. The dashed line 
corresponds the Coulomb interaction. The wave line corresponds the fine or hyperfine interaction. 
G is the reduced Coulomb Green's function. 



This equation contains the relativistic correction of order {Za)^ , which can be calculated 
on the basis of the Dirac equation 0, H^, the correction of order ai^Za)^ enhanced by the 
In(Za) BHIIJ, a number of terms of fifth and sixth order over a which are determined 
by the effects of the vacuum polarization |^. The relativistic recoil effects of order 
mi iZa)^ Irri'). in the energy spectra of hydrogenic atoms were investigated in Refs. [l|, 0, 
|39|, |40|. In the fine splitting (3) they were calculated in [34, 40]. Additional corrections of 
the same order were obtained in |4lJ. They don't depend on the muon total momentum j 
and give the contribution only to the Lamb shift. The contributions to the coefficients Ayp 
and Byp are specified by the first and second order perturbation theory. Numerical values 
of the terms in the expression (3) presented in the analytical form, are depicted in Table I 
with the accuracy 0.001 ^eV . Such precision is related to the numerical values of a number 
of contributions to the energy spectrum obtained below from the potentials containing VP 
effects. The fine structure interval (3) in the energy spectrum of electronic hydrogen is 
considered for a long time as a basic test of quantum electrodynamics [i^, 4^, 3 • 

The interaction operator which gives the contribution to the coefficient Ayp-, is repre- 
sented by the Feynman diagrams in Fig.l. The vacuum polarization effect leads to the 
modification both the Coulomb interaction and the spin-orbit interaction in expressions (1), 
(2): 
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l(2s^+l)/s^, me is the electron mass. The diagram 
(a) in Fig. 1 gives the first order perturbation theory contribution to the fine structure of 
order a{Za)^. Averaging the potential (2) over the wave functions of the 2P - state 
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we obtain the following contribution to the interval (3): 
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Notice that the higher order corrections a'^{Za)'^ entering the are taken into account 
in this expression as well as the recoil effects. The same order contribution a{Za)^ can 
be obtained in the second order perturbation theory (see the diagram (b) in Fig.l). In 
the second order perturbation theory the energy spectrum is determined by the reduced 



Coulomb Green's function whose partial expansion has the form |44 
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The radial function can be represented in the form of the Sturm expansion in the Laguerre 
polynomials. For the 2P- state this function is the following: 
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+30^(4440 - 4230^' + 1152/^ - 109^'^ + 3^'^ 
where z = fiZar, are the Laguerre polynomials defined by the relation: 
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To use the expression (9) in specific calculations [27|, |29|, |30| we have to integrate over z, 
z' and after that to sum over m. Another representation for the G2p{r,r') was obtained in 
Ref.jisjl by summing over m in (9): 
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where z^ = min{z, z'), z^ = max{z, z'), C = 0.577216... is the Euler constant, z = Wr. 
Using Eqs. (11) and (12) we transform the second part of the correction (the first part 
is given by Eq.(7)) of order a{ZaY to the fine structure appearing in the second order 
perturbation theory as follows: 
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Let us consider the two-loop vacuum polarization contributions in the one-photon inter- 
action shown in Fig. 2. They give the corrections to the fine splitting of P- levels of order 
a\Za)\ 



a b c 

FIG. 2: Effects of two-loop electron vacuum polarization in the one-photon interaction. 



To find the interaction operator in the momentum representation for the amphtude (a) 
in Fig. 2 it is necessary to make the double change 
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in the photon propagator. Then in the coordinate representation the interaction operator 
takes the form: 
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Averaging (15) over wave functions (6), we obtain the correction to the interval (3): 
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The two-loop vacuum polarization amplitudes in the diagrams (b,c) of Fig. 2 can be calcu- 
lated using the modification of the photon propagator of the different type [45l |: 
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To obtain the numerical value of the two-loop contribution in this case it is convenient to 
use the coordinate representation in which the potential corresponding to the amplitudes 
(b,c) in Fig.2 has the form: 
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Its contribution to the fine splitting (2P3/2 — 2P1/2) in muonic hydrogen can be written in 
the integral form: 
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FIG. 3: Effects of two-loop electron vacuum polarization in the second order perturbation theory. 
The dashed line corresponds to the Coulomb interaction. The wave line corresponds to the fine or 
hyperfine interaction. G is the reduced Coulomb Green's function. 

Two-loop vacuum polarization contributions in the second order perturbation theory 
shown in Fig. 3, have the same order Q^{ZaY. For their calculation we can use earlier 



obtained relations (2), (4), (5 ), ( ll)] and the modified Coulomb potential because of the 

iz; 

2 roo roo /' Za'- 



two- loop vacuum polarization |29l I30l|: 
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The amplitude (a) in Fig.3 gives the following correction of order a'^{Za)^ to the fine split- 
ting: 
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Two other contributions from the amphtudes (b), (c) in Fig.3 have the similar integral 
structure. Their numerical values are included in Table I. The summary result for the fine 
splitting AE^'^ in (/ip) is presented here also. It accounts the numerous earlier performed 
calculations discussed in the review article [l| and new corrections obtained in this work. 



TABLE I: Fine structure of P-wave energy levels in muonic hydrogen. 



Contribution to fine 
splitting AS^^ 


Numerical value 
of the contribution in fieV 


Reference, 
equation 


Contribution of order {Za)'^ 

t?{Zat , mi \ 
32m^ V 277*2^ 


8329.150 


[21, 23], (3) 


Muon AMM contribution 

i6m^ V ; 


17.637 


[21, 23], (3) 


Contribution of order {Za)^: ^™25^""* 


0.312 




Contribution oi order [Za) mi/m2- Q4my 


-0.028 


(3), [34, 40] 


Contribution of order a{Za)^ 

• J 1 J 1 T ~\ I I \ /ATT' \ 

m the nrst order Fi {AVyp} 


3.046 


(7), [21, 23] 


Contribution of order a{Za)^ 
m the second order Fi 
{AVyp ■ G ■ AVf') 


1.928 


(13) 


Contribution of order a{Za)^ 

a(Za)''u^ Ft u(Za)^ . ll 
32nm'i [ mi ~'~ 5j 


-0.008 


fi] 


VP Contribution in the second 
order FT of order a'^{Za)'^ 
(AV^p-G-AV^'p) 


0.002 


(23) 


VF Contribution from I'j interaction 
of order a'^{Za)^ {AV^p_yp) 


0.002 


(16) 


VF Contribution from I7 interaction 
of order a^Za)^ {^vtioopyp) 


0.021 


(20) 


VF Contribution in the second 
order FT of order a^{ZaY 
{AV^P^yp-G-AVf^) 


-0.002 


(21) 


VF Contribution in the second 
order FT of order a'^{ZaY 
{^y^-ioovyp ■ G ■ AVf^) 


0.022 


(22) 


Summary contribution 


8352.082 





III. HYPERFINE STRUCTURE OF THE ENERGY LEVELS 2P1/2 2P3/2 

The leading order contribution to the hyperfine splitting of the energy levels 2Pi/2 and 
2P3/2 in muonic hydrogen of order {Za)^ is determined by the following potential (the 
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hyperfine part of the Breit potential) [33|: 
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where n = r/r, k is the proton anomalous magnetic moment. The operator (24) does not 
commute with the operator of the muon total angular momentum J = L + \cri- So, the 
energy levels 2Pi/2 and 2P3/2 are mixed and the hyperfine structure of P-wave levels is more 
complicated. 

To calculate the diagonal matrix elements (2P1/2I AV^''^'*|2Pi/2) and (2P3/2I AV^-^*|2P3/2) 
we can use the following replacement for the operators (S1S2) (LS2) containing the nuclear 
spin: 
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where (siJ), (LJ) are eigenvalues of corresponding operators between the states with equal 
orbital momentum /. Moreover, the angle averaging in the second term of (24) can be carried 
out by means of the relation |46[: 



{5ij - MiUj) = ~{4:6ij - 3LiLj - SLjLi). 
The diagonal matrix elements have the following general structure: 
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where Ep = {Za)'^n^{l + /t)/3mim2 is the Fermi energy for the level with the n = 2. The 

1/2 3/2 

calculation of the relativistic corrections A,/^i , A^^i in this formalism includes the study 
of two-photon, three-photon exchange diagrams and the second order perturbation theory 
contributions with the Breit Hamiltonian (1), (2), (24). More simple approach to their 
calculation is based on the Dirac relativistic equation j37|, |47|. In this case the hyperfine 
splitting potential has the form: 
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and gives the following contributions to the hyperfine structure: 
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where the nuclear magnetic moment /j, = gNfJ'N^2 if^N = e/2m2). The typical radial integrals 
Rk = Io° dkfkdr for this case are determined by the Dirac wave functions fk,gk of the 
states 2Pi/2, 2P3/2. Accounting their exact form |46| we obtain the following values of the 
relativistic corrections to the hyperfine structure of the P - wave level: 
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These values of the coefficients coincide with the analytical calculation of the contribution 
of order m^(Za)^/m2 in the hyperfine structure of the P - levels of hydrogen atom for the 
n = 2, carried out in Ref . [i^ . 

The fifth order contribution over a appears in the hyperfine splitting as well as in the 
fine structure due to the electron vacuum polarization (see diagrams (a), (b) in Fig.l). The 
modification of the hyperfine part of the Breit potential from the vacuum polarization is 
described by the following relation (the substitution (14) is used) (23l |: 
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For the subsequent transformations of the diagonal matrix elements of the operator (32) we 
use the relation for the angle average (26) and the similar expression for the first term in 
the square brackets of Eq.(32): 



(Sij - riiTij) = i[25y + LiLj + LjLi]. 



(33) 



Then the contributions of the vacuum polarization in the first and second orders PT can be 
written as follows: 
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It should be noted that the sums of the corrections (34), (36) and (35), (37) in the terms 
Aypa, Aypa are equal 0.00025 and 0.000056 correspondingly. They are slightly different 
from the results 0.00022 and 0.00008 obtained in [23^1. Neglecting the recoil effects and muon 
anomalous magnetic moment in the expressions (34) and (35) we obtain the results which 
coincide with the calculation based on the Dirac equation with the potential: 

A^y^D(^) = e^^^^T^ rp{s)ds{l + 2m,sr)e-^^^'\ (38) 

\l1 3/2 

Numerical values of the contributions AypOt and Aypa in this case are equal 0.00019 and 
0.00002 correspondingly. 



TABLE II: Hyperfine structure of P-wave energy levels in muonic hydrogen. 



Contribution to hyperfine 


Numerical value 


Numerical value 


Reference, 


splitting 


01 the contribution m 
A^«/^(2Pi/2), iieV 


01 the contribution m 
A^'^/^(2P3/2), neV 


equation 


Contribution of order {Zoif" 


7953.195 


3392.112 


(27), (28), [21, 22] 


Muon AMM contribution 


4.432 


-0.886 


(27), (28), [21, 23] 


Contribution of order (Zoif' 


1.092 


0.065 


(31)J401 


Contribution of order a(Za)'* 


3.830 


0.497 


(34),(35) 


in the first order PT [ISVy^p) 






[21, 2^ 


Contribution of order a(Za)'* 






(36),(37) 


in the second order PT 


1.838 


0.783 




(AVVp • G ■ AV^f') 








VP Contribution in the second 






(4),(32) 


order PT of order a'^{Za)'^ 


0.002 


0.002 




{AV^P-G-AV^{.') 








VP Contribution of I7 interaction 


0.003 


0.00004 


(39) 


of order a'^{Za)^ {AVy^p_yp) 








VP Contribution of I7 interaction 


0.026 


0.008 


(40) 


of order a2(Za)4 {AV^H^ yp) 








VP Contribution in the second 


-0.001 


-0.0003 


(11),(21),(24) 


order PT of order a^{Za)^ 








{AV^p.yp ■ G ■ AV^f') 








VP Contribution in the second 


0.017 


0.007 


(11),(22),(24) 


order PT of order a^(Za)^ 








{^V,^-ioo,yp-G-AV^f^) 








Summary contribution 


7964.364 


3392.588 





Two-loop vacuum polarization corrections in the hyperfine part of the potential for the 
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1.278 fieV. 



11 



/ 7^ obtained by means of Eqs.(14), (17) have the integral representation: 
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Omitting further details of the calculation of their expectation values which can be performed 
analogously (20) and (23), we represent in Table II numerical results of the contributions 
of the potentials (39), (40). Another part of two-loop corrections of order oP'iZoi)'^ to the 
hyperfine structure in the second order PT is shown in Fig. 3. We also included numerical 
results obtained by means of these amplitudes in Table II. 

Nondiagonal matrix element has an important role to attain the high accuracy of the 
calculation of the P-wave levels in muonic hydrogen. We present its general structure as 
follows: 
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48 , 



mi(l + 2/«) 
m2(l + K) 



P" 



(41) 

where for the simplicity we have restricted by the fifth order terms over a in the vacuum 
polarization effects, the fifth and higher order terms in the muon anomalous magnetic mo- 
ment, the relativistic effects of order [Zoif and the recoil effects. First three terms in the 
right part of Eq. (41) appear from the potential (24). In the Dirac theory the relativistic 
corrections are determined by nondiagonal radial integral: 



Pi 3 

2 2 



(c/i(r)/|(r) + ^|(r)/i(r))c/r. 



(42) 



After the integration inj_42) which can be performed using explicit form of the wave functions 
/i/2,3/2('^)5 5'i/2,3/2('^) |46|, we obtain the coefficient C^ei = 9/16. To calculate the vacuum 
polarization effects we substitute the potential (32) into 7. Then we have to calculate the 
matrix elements of the following operators: 



Ti = (Lo-2), T2 = [cr^cr^ - 3(o-^n) (0-211)] , T3 = [cr^a^ - {a^n)(a^n)] 



(43) 
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FIG. 4: The structure of S'-wave and P-wave energy levels in muonic hydrogen for the n = 2. 

After the angle averaging in Eq.(43) by means of Eqs.(26) and (33) these matrix elements 
can be expressed in terms of the 6j - symbols: 



where the value of the total momentum F = 1 (F = S2 + J), / = 1,_L,= \/2j~+T, / = 
A/2j' + 1, and the numerical values of the 6j - symbols are taken from |4S|. As a result the 
vacuum polarization contributions to the nondiagonal matrix element (41) in the first and 
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second order PT have the form: 



71 



2T7ie s \ 

w ) X 



X 



mi(l + 2k) 
2m2(l + k) 



72 = f Pi/2|Alfp • G ■ Ay^^YP3/2) = Ef 



2meS Arn^s^ 
-X tttt;— a; 



V2 



W 



a 



X 



p{s)ds 



'tic d/Oc 6 



-X i+- 



2592 / TT 



1 + 



mi(l + 2k 
^2(1 + k) 



(45) 

0.617 peV, 
X (46) 



X' 



/2 



-0.184 /ieV. 



Putting = 0, mi/m2 = in Eq.(45), we obtain the value —0.889 fieV which coincides 
with the calculation of this matrix element in the Dirac theory with the potential (38). 
Summary numerical value of the matrix element (41) is equal 7 = —796.192 fieV. It leads 
to the shift of the energy levels 2'^P3/2 and 2^Pi/2 by the value 6 = 144.560 fieV as shown 
in Fig. 4. 



IV. SUMMARY AND CONCLUSION 

In this work we calculate the QED effects in the fine and hyperfine structure of the 2Pi/2, 
2P3/2 energy levels in muonic hydrogen. The electron vacuum polarization contributions of 
orders a^, with the recoil corrections, the relativistic effects of order are considered. 
Numerical values of the contributions are listed in Tables I and II. We give the references 
on the papers devoted to the leading order calculation of the structure of P-wave levels. 
For the comparison of the obtained results with the earlier performed calculations we used 
the review article containing the modern status of the investigations in the physics of 
relativistic energy spectra of simple atoms. 

Let us summarize the basic particularities of the calculation performed above. 

1. Special attention in our investigation has been concentrated on the vacuum polar- 
ization effects. For this purpose we obtain the terms of the interaction operator in muonic 
hydrogen which contain the one-loop and two-loop vacuum polarization corrections. 

2. At each order over a we retain the recoil effects in the terms ~ mi/m2- The experi- 
mental values of the muon and proton anomalous magnetic moments are used jif. 

3. The calculation of the relativistic corrections to the diagonal and nondiagonal matrix 
elements both for the fine and hyperfine structure intervals is performed by means of the 
Dirac equation. In the second order perturbation theory we use the compact representation 
for the reduced Coulomb Green's function obtained in Ref . jisl] . 

Total numerical values for the fine structure interval AE^^ (3) and hyperfine structure 
splittings of 2P1/2, 2P3/2 states are presented in Tables I, II. Accounting also our calculation 
of the mixing between states 2^P3/2 and 2^Pi/2 (the correction 5), we obtain the change 
of the hyperfine splittings on the 6 = 144.560 j^eV: AE^f''{2Pi/2) = AE^f'{2Pi/2) - 6 = 
7819.804 fieV, AE^f'{2P^i2) = AE^f'{2P3/2) - 6 = 3248.028 fieV. The theoretical error of 
the obtained results is determined by the contributions of higher order and amounts up to 
10^^. The results of this work improve the previous calculations in (2ll . lisl . l26l | because of 
the investigation the effects of order and can be considered as a reliable estimate for the 
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fine and hyperfine structure intervals for tfie P- levels in muonic hydrogen. These results are 
important for the experiment at PSI jisl . [l9| . The disposition of the P-wave energy levels, 
shown in Fig. 4, is determined by the following numerical values: —5973.27 ^eV (2^Pi/2), 
1846.53 /zel/ {2^Pi/2), 6376.28 ^leV {2^Pz/2), 9624.30 ueV" (2^^3/2)- Accounting the Lamb 
shift value {2P — 25*) in (yup), obtained in Refs.0, Ho], the hyperfine sphtting of the 25*- 
state from [2^ and the results of present study we obtain the value of the energy interval in 
muonic hydrogen ^E{2^P^/2 - 2^81/2) = 205975.6 fieV. 
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